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Summary 

Our  team  proposed  the  development  of  a 
modular,  adaptive,  parallel  software  framework 
for  modeling  the  Sun-Earth  system.  The  goal  of 
the  project  was  the  development  of  a  large-scale 
model  of  the  solar-terrestrial  environment 
allowing  a  fuller  understanding  of  space  weather 
and  a  framework  to  test  theories  and  investigate 
the  broad  implications  of  new  observations. 
Particular  attention  was  to  be  devoted  to  CME 
formation,  propagation,  and  interaction  with  the 
magnetosphere;  SEP  acceleration  in  the  low 
corona,  SEP  acceleration  in  the  interplanetary 
medium,  and  SEP  transport. 

The  proposed  model  was  to  be  driven  by  real¬ 
time  observations  that  are  incorporated  into  a 
validated  set  of  tightly-coupled  models.  We 
proposed  that  the  model  would  be  developed  as  a 
set  of  modular  routines  capturing  the  physics  of 
interacting  domains  in  the  space  environment, 
providing  flexibility  for  revision  and  improvement 
of  the  modules.  The  proposed  new  space  weather 
modeling  framework  was  envisioned  to  serve  four 
complementary  goals: 

1,  Provide  a  flexible  "plug-and-play'  software 
framework  for  DoD’s  space  weather  research. 
Models  of  physics  domains  and/or  processes 
can  be  "plugged"  into  the  framework  and 
researchers  can  lest  the  response  of  the  entire 
Sun-Earth  system  to  the  new  modules. 

2.  Enable  the  proposing  team  to  develop,  test 
and  validate  improved  models  of  solar 
explosive  event  initiation  (such  as  CMEs), 


interplanetary  propagation,  SEP  generation 
and  gcoeffectiveness. 

3.  Provide  a  complete  Sun-to-Earth  space 
weather  simulation  tool.  When  run  from  real¬ 
time  observations  of  the  region  near  the  Sun, 
the  solar  and  interplanetary  modules  would 
have  the  capability  to  provide  an  up-to-date 
display  of  the  calculated  present  state  of  the 
solar  wind  in  the  inner  heliosphere,  as  well  as 
a  prediction  for  the  next  24-48  hours  near 
Earth. 

4.  The  model  would  be  made  available  to  the  Air 
Force  and  NOAA  SEC  as  a  prototype  for 
validation  and  operational  testing. 

We  are  very  proud  to  note  that  all  four  goals 
have  been  met.  In  addition,  this  project  supported 
a  large  number  of  science  publications, 
conference  presentations  as  well  as  Ph.D.  students 
and  postdocs.  In  the  body  of  this  final  report  we 
provide  details  of  the  Space  Weather  Modeling 
Framework,  some  of  the  science  highlights  and 
programmatic  experience  we  gained  with  this 
exciting  project. 

The  COSTEM  project  also  supported  over  ten 
Ph.D.  students  during  its  6  years  of  funding.  It 
also  helped  3  young  scientists  to  achieve  national 
prominence  and  obtain  highly  coveted 
instructional  faculty  positions  (ilia  Rousscv. 
Aaron  Ridley  and  Mike  Liemohn). 

Overall,  COSTEM  was  a  successful  project. 
We  all  wish  we  could  have  more  projects  like  this. 
We  published  over  100  peer  reviewed  papers, 
gave  more  than  100  invited  and  174  contributed 
presentations  at  professional  meetings. 
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1  Introduction 

The  PI  has  been  active  in  space  science  for 
over  three  decades*  but  this  MUR  I  grant  was  one 
of  the  most  enjoyable*  productive  and  exciting 
projects  he  has  ever  participated  in.  We  would 
like  to  express  our  thanks  to  the  DoD  for  giving  us 
the  opportunity  to  cany  out  this  research,  to 
AFOSR  for  managing  the  project,  and  to  USAF 
Majors  Dr.  Paul  Bcllaire  and  Dr.  David  Byers  for 
providing  a  supportive  and  positive  environment 
and  management* 

The  main  goal  of  the  COST  EM  project  was  to 
create  a  uscr-fricndly  modeling  tool  that  enables 
the  DoD  and  scientific  communities  to  simulate* 
and  eventually  forecast,  the  system  behavior  of 
the  Sun- Earth  space  environment. 

A  very  important  factor  in  the  success  of  the 
CQSTEM  project  was  the  synergism  with  our 
NASA  funded  Computational  Grand  Challenge 
(CGC)cfforl  that  supported  the  development  of 
the  computational  infrastructure  of  the  Space 
Weather  Modeling  Framework  (SWMF)  in  close 
collaboration  with  the  Earth  System  Modeling 
Framework  (ESMF).  The  synergism  between 
COSTEM,  SWMF*  and  ESMF  was  additionally 
enhanced  by  an  NSF  Information  Technology 
Research  (ITR)  grant  that  was  providing  support 
for  research  in  data  assimilation  methodologies 
applicable  to  the  space  environment,  as  well  as 
^threat  adaptive”  grid  computing  technologies, 
where  we  dynamically  allocate  distributed 
computational  resources  and  increase  the 
resolution  and  throughput  of  the  simulation  on  the 
fly  when  a  threatening  space  storm  occurs  on  the 
Sun. 

The  synergism  between  our  MUR1,  CGC  and 
ITR  projects  was  made  possible  by  the  support  of 
the  science  managers  at  NASA,  NSF  and  AFOSR: 
Rich  Behnke,  Paul  Bellairc,  Joe  Bredekamp, 
David  Byers  and  Jim  Fischer  Wc  made  their  job 
easier  by  clearly  identifying  and  separating  the 
three  projects  at  Michigan:  NASA  funds  were 
used  to  develop  computational  technology,  NSF 
funds  were  used  for  data  assimilation*  grid 
computing  and  model  validation,  while  DoD 
funds  were  devoted  to  model  development  for  the 
3D  global  corona*  CME  initiation,  solar  energetic 
particles*  inner  magnetosphere,  ionosphere- 
thermosphere*  radiation  belts,  and  other  regional 
phenomena. 


The  integration  of  the  three  large  projects 
represented  considerable  management  challenges 
(we  will  discuss  these  issues  later),  but  it  also 
gave  us  a  unique  opportunity  to  create  something 
that  is  more  than  the  sum  of  its  parts.  The 
resulting  model  suite  is  the  most  advanced 
computational  tool  available  today  to  simulate  the 
complex  behavior  of  the  Sun-Earth  space 
environment.  It  is  a  unique  tool  that  has  been 
extensively  used  by  our  team  for  scientific 
research  (wc  published  over  100  peer  reviewed 
papers  over  the  lifetime  of  our  MUR  I  project).  In 
addition,  it  is  continuously  running  at  CCMC  in 
real  time  as  a  precursor  to  operation  transition  to 
NO  A  A  and  the  Air  Force. 


In  the  following  sections  we  discuss  the  main 
accomplishments  of  the  COSTEM  project,  we 
highlight  some  science  results*  outline  new 
investigations  that  were  made  possible  by  the 
success  of  the  COSTEM  project,  describe  the 
management  strategy,  and  most  importantly,  wc 
discuss  the  lessons  learned  from  this  exciting 
undertaking. 


Page  2  of  32 


MUR!  F4 9 620*01-1 -0359  Final  Report 


PI:  Tam  as  Gombosi 


2  ACCOMPLISHMENTS 

2*1  The  Four  Goals  of  COSTEM 

2.1.1  "Plug-and-Play”  Software  Framework 
SWMF 

The  Space  Wcaihcr  Modeling  Framework 
(SWMF)  [29,  43]  was  designed  in  2001  and  has 
been  developed  to  integrate  and  couple  several 
independently  developed  space  physics  codes  into 
a  flexible  and  efficient  model  Figure  I  shows  the 
schematics  of  the  SWMF  and  its  components.  The 
main  design  goals  of  the  SWMF  were  to  minimize 
changes  in  the  original  models,  provide  a  general 
flexible  and  efficient  method  to  execute  the 
coupled  models  on  massively  parallel  distributed 
memory  machines,  and  to  allow  adding  new 
components  and  new  physics  models  for  existing 
components  with  ease  [29], 

The  SWMF  is  written  in  Fortran  90  and  Perl 
and  has  been  designed  in  an  object-oriented 


manner  from  the  beginning.  Data  naming 
conventions  and  the  use  of  CVS  were  adopted 
from  the  beginning.  Rigorous  testing  procedures 
have  been  applied  from  the  beginning  as  well. 
Most  of  the  implemented  modules  contain  unit 
tests.  The  components  are  tested  individually,  and 
the  functionality  of  the  whole  SWMF  is  tested  by 
system  tests  exercising  multiple  coupled 
components.  The  SWMF  and  its  components  pass 
this  comprehensive  hierarchical  test  suite  every 
night  on  several  computcr/eompilcr  platforms. 

The  SWMF  source  code  is  well  documented. 
The  reference  manual  is  generated  from  the 
source  code  using  Perl  scripts.  All  input 
parameters  arc  documented  in  XML  files,  which 
are  converted  into  the  user  manual  We  have  also 
developed  a  graphical  user  interface  (SWMF 
GUI)  to  improve  the  usability  of  this  complex 
software. 

We  are  regularly  extending  the  SWMF  with 


http://cs6vn.0ngl  n.umlch  .sd  u 

Figure  I  Schematic  of  coupling  in  SWMF 
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new  components  and  new  component  versions.  It 
takes  typically  two  weeks  of  work  from  two 
developers  to  add  a  new  component  to  the 
SWMF.  Most  of  this  time  is  spent  on  adopting  the 
physics  model  to  the  SWMF  and  on  writing  and 
testing  the  couplers  between  the  components. 
New  components  that  do  not  use  Cartesian  grids 
can  take  longer  if  they  require  high-performance 
grid  rendezvous  algorithms.  The  SWMF  currently 
covers  10  physics  domains,  many  of  them 
represented  with  multiple  models.  The  SWMF  is 
capable  of  simulating  the  Sun-Earth  system  from 
the  solar  corona  to  the  upper  atmosphere  of  the 
Earth  faster  than  real  time  on  hundreds  of 
processors. 

BATS-R-US 

The  Block  Adaptive  Tree  Solar-wind  Roe 
Upwind  Scheme  [cf.  60]  has  been  designed  and 
developed  to  solve  magnetohydrodynamic  (MHD) 
space  physics  problems,  BAT5-R-US  was 
intended  from  the  beginning  to  use  massivdy- 
parallel,  distributed-memory  machines  efficiently. 
Therefore,  it  was  decided  that  the  code  would  be 
written  from  scratch  using  Fortran  90,  the  MPI 
library  for  communication,  and  a  block -adaptive, 
Cartesian  grid  structure  [60].  The  resulting  code 
showed  excellent  efficiency,  scaling  properties 
(see  Figure  2)  and  portability,  and  resulted  in 
state-of-the-art  space  science  simulations. 

In  the  past  12  years,  BATS-R-US  has  been 
continuously  developed  in  several  ways.  New 
numerical  algorithms  were  added:  Rusanov,  HLL, 
and  HLLD  solvers,  and  various  techniques  to 
control  the  divergence  of  the  magnetic  field  [92], 
implicit  time  stepping  [23]  parallel  field  line 
tracing  [52],  oscillation-free  second-order 
schemes  at  resolution  changes  [16],  etc.  The 
physical  equations  were  extended  from  ideal 
MHD  to  resistive  MHD,  semi-relativistic  MHD 
[91  ],  multi-species  MHD  and  very  recently  to  Hall 
MHD  and  multi-fluid  MHD.  We  note  that  the 
MHD  solvers  can  be  used  for  the  equations  of 
compressible  hydrodynamics  without  any 
modification,  The  block -adaptive  grid  can  now 
use  generalized  coordinates,  including  spherical, 
cylindrical  and  toroidal  grids. 

With  better  and  more  flexible  algorithms,  the 
range  of  applications  has  also  widened 


dramatically:  BATS-R-US  has  been  used  to  model 
almost  all  planets,  several  moons,  and  comets  of 
the  solar  system,  flux  emergence  from  the 
photosphere  to  the  corona  [64],  coronal  mass 
ejections  [63,  75  j,  the  solar  corona  [70],  the  outer 
heliosphere,  etc.  These  applications  imply  that 
BATS-R-US  has  been  extensively  validated 
against  observational  data.  The  code  also  enables 
one  to  use  special  grids,  such  as  spherical,  toroidal 
or  the  vertex-based  stretched  grids. 

As  the  complexity  of  BATS-R-US  increased, 
more  and  more  of  the  standard  software 
engineering  practices  were  adopted.  Wc  have 
been  using  the  concurrent  version  control  system 
(CVS)  for  several  years.  Much  of  the  code  is  now 
written  in  a  modular  and  object  oriented  fashion. 
The  developers  introduced  a  data  naming 
convention  so  the  code  is  more  readable  and 
easier  to  maintain.  The  functionality  of  BATS-R- 
US  is  tested  with  a  suite  of  comprehensive  tests 
run  nightly  on  half  a  dozen  different  computers 
and  compilers.  The  input  parameters  arc 
documented  and  described  in  XML  format  and 
checked  with  a  script  for  correctness.  Most  of  the 
documentation  is  generated  automatically. 

2.1.2  Development  Testbed 

With  the  help  of  the  SWMF-bascd  testbed  wc 
were  able  to  develop  and  apply  several  new 
models  for  the  quiet  solar  wind  (synoptic  wind), 
CME  initiation  mechanisms,  heliospheric 
transient  propagation,  SEP  models,  and  the 
interaction  of  interplanetary  transients  with  the 
magnetosphere-ionosphere  system. 
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Synoptic  Solar  wind, 

A  prerequisite  of  any  successful  end-to-end 
space  weather  model  is  the  ability  to  simulate  the 
background  solar  wind  filling  the  inner 
heliosphere  at  the  time  of  the  solar  eruption.  The 
ambient  conditions  are  crucial  for  the  successful 
modeling  the  propagation  of  a  time  dependent 
phenomena  (such  as  CMEs), 

We  developed  a  new  MHD  model  for 
simulating  the  large-scale  structure  of  the  solar 
corona  and  solar  wind  under  “steady  state*’ 
conditions  stemming  from  the  Wang-Sheclcy- 
Arge  empirical  model  and  driven  by  high- 
resolution  MDI  solar  magnetograms  [5,  70],  The 
processes  of  turbulent  healing  in  the  solar  wind 
were  parameterized  using  a  phenomenological, 
thermodynamical  model  with  a  variable 
polytropic  index.  We  employed  the  Bernoulli 
integral  to  bridge  the  asymptotic  solar  wind  speed 
with  the  distribution  of  the  polytropic  index  on  the 
solar  surface.  This  “synoptic  wind"  model 
successfully  reproduces  the  bulk  properties  of  the 
“quiet"  solar  wind  at  Earth  for  solar  minimum 
conditions  (see  Figure  3). 

This  is  significant  progress  over  previous 
work,  but  obviously  more  work  is  needed  to  have 
a  validated,  reliable  synoptic  wind  model  that  can 
be  used  for  all  solar  cycle  conditions. 

CME  INITIATION 

One  of  the  very  active  research  areas  in 
hcliophysics  is  CME  initiation.  This  is  a  very 
important  unsolved  question  that  has  profound 
consequences  for  space  weather  events.  With  the 
help  of  the  SWMF-bascd  simulation  testbed  we 
investigated  several  magnetically  driven  CME 
initiation  mechanisms. 

We  modeled  the  emergence  of  a  magnetic 
flux  rope  passing  from  below  the  photosphere  into 
the  corona  [59],  For  the  initial  state,  wc  prescribed 
a  plane-parallel  atmosphere  that  comprises  a 
poly  tropic  convection  zone,  photosphere, 
transition  region,  and  corona.  Embedded  in  this 
system  was  an  isolated  horizontal  magnetic  flux 
rope  located  10  photospheric  pressure  scale 
heights  below  the  photosphere.  The  flux  rope  was 
uniformly  twisted,  with  the  plasma  temperature 
inside  the  rope  reduced  to  compensate  for  the 
magnetic  pressure.  Density  was  reduced  in  the 
middle  of  the  rope,  so  that  this  section  buoyantly 


rose.  The  early  evolution  proceeded  with  the 
middle  of  the  rope  rising  to  the  photosphere  and 
expanding  into  the  corona.  Just  as  it  seemed  the 
system  might  approach  equilibrium,  the  upper  part 
of  the  flux  rope  began  to  separate  from  the  lower, 
mass-laden  part.  The  separation  occurred  through 
stretching  of  the  field,  which  forms  a  current 
sheet,  where  reconnection  severs  the  field  lines  to 
form  a  new  system  of  closed  flux.  This  flux  then 
erupts  into  the  corona.  Essential  to  the  eruption 
process  were  shearing  motions  driven  by  the 
Lorcntz  force,  which  naturally  occur  as  the  rope 
expands  in  the  pressure-stratified  atmosphere.  The 
shearing  motions  transport  axial  flux  and  energy 
to  the  expanding  portion  of  the  magnetic  field, 
driving  the  eruption. 

We  also  developed  another  CME  initiation 
model  [75]  that  was  considering  the  loss  of 
equilibrium  of  the  three-dimensional  flux  rope 
configuration  of  Titov  &  Demoulin,  We  arc  able 
to  determine  the  conditions  for  which  stable 
equilibria  no  longer  exist.  Our  results  imply  that  it 
is  possible  to  achieve  a  loss  of  equilibrium  even 
though  the  ends  of  the  fiux  rope  are  anchored  to 
the  solar  surface.  However,  in  order  to  have  the 
fiux  rope  escape,  it  is  necessary  to  modify  the 
configuration  by  eliminating  the  arcade  field. 


SEP  Transport 

Solar  energetic  particles  primarily  propagate 
along  interplanetary  magnetic  field  lines 
connecting  the  solar  corona  to  the  interplanetary 
medium.  Existing  SEP  acceleration  and  transport 
models  use  prescribed  magnetic  field  lines  to 
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investigate  the  transport  and  diffusive  shock 
acceleration  of  solar  energetic  particles. 

With  the  help  of  the  SWMF  testbed  we  were 
able  to  carry  out  the  first  self-consistent 
simulation  of  SEP  events  by  coupling  field- 
aligned  SEP  transport  models  to  simulated 
magnetic  field  lines  extracted  from  a  time- 
evolving  simulation  of  a  CME  transient 
propagating  through  the  heliosphere.  The 
technique  is  based  on  a  new  field  line  advection 
model  [57]  that  is  now  part  of  SWMF.  The 
transport  of  SEPs  can  now  be  simulated  by 
different  SEP  transport  models  [25,  45,  61]  thus 
enabling  us  to  look  for  experimentally  verifiable 
distinguishing  predictions. 

2.1.3  Space  Weather  Simulation  Tool 

The  SWMF  with  its  existing  model  suite  is  an 
end-to-end  space  weather  simulation  tool.  It  has 
been  used  to  model  space  weather  events  from  the 
Sun  to  Earth  [1,  19,  33,  41,  63,  100].  These 
simulations  employed  increasingly  sophisticated 
models  and  approaches  to  simulate  the  Sun-Earth 
system.  In  our  latest  simulation  [1]  we 
successfully  reproduced  many  signatures  of  the 
most  geoeffective  Halloween  storm  (October  29; 
2003).  In  particular,  we  not  only  succeeded  in 
reproducing  the  main  features  of  the  CME  at 
Earth,  but  were  also  successful  in  matching 
magnetospheric  observations  in  the  dayside 
magnetosphere,  in  the  nightside  closed  field  line 
region,  and  in  the  distant  magnetotail. 


Figures  5  shows  a  comparison  between 
simulated  and  observed  white  light  images  shortly 
after  the  CME  eruption.  We  note  the  excellent 
agreement  between  the  shapes  of  the  observed  and 
simulated  structures,  as  well  as  the  location  of  the 
shock.  This  is  a  very  encouraging  progress  in 
modeling  large  CMEs  as  they  propagate  from  the 
low  corona  to  Earth. 

Figure  6  compares  observed  and  simulated 
plasma  parameters  in  the  magnetosphere.  We  note 
the  excellent  overall  agreement.  It  is  particularly 
challenging  to  match  the  magnetospheric 
parameters  simultaneously  in  distinct  regions  in 
the  magnetosphere.  At  the  time  of  this  event  the 
Cluster  spacecraft  was  in  the  dayside 
magnetosphere.  Polar  was  magnetically  connected 
to  the  high  latitude  region,  GOES- 10  was  behind 
the  Earth  in  the  closed  field  line  region,  and  Wind 
in  the  distant  magnetotail  some  150  R[ 
downstream. 


2. 1.4  Transitioning 

One  of  the  primary  goals  of  our  MURI  project 
was  to  create  a  space  weather  simulation  tool  that 
is  robust,  easy  to  use  (well,  relatively)  and  can  be 
eventually  transitioned  to  operations.  An  early 
attempt  to  transition  BATS-R-US  to  Air  Force 
operations  was  unsuccessful,  because  the  lack  of 
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trained  personnel  at  the  operation  center  the 
complicated  user  interface  and  the  lack  of  daily 
communication  between  the  code  developers  and 
the  operators. 

We  also  had  discussions  with  NOAA  SEC 
about  transitioning  SWMF,  The  new  director  of 
SEC,  Dr,  Thomas  Bogdan  is  extremely  interested 
in  running  SWMF  in  real  time  and  validating  the 
results.  However  presently  SEC  does  not  have  the 
computational  resources  to  be  able  to  run  SWM  in 
real  time  (or  even  in  near  real-time).  In  the  next 
year  or  so  they  will  have  access  to  larger  NOAA 
computing  resources  and  SEC  will  transition  the 
latest  version  of  SWMF  and  start  real-time  testing. 

The  interagency  Community  Coordinated 
Modeling  Center  (CCMC)  was  the  First  external 
user  of  BATS-R-US  and  SWMF.  CCMC  has  the 
personnel,  the  expertise  and  the  computing 
resources  to  utilize  our  high  performance  code. 
Wc  are  in  regular  communication  with  CCMC 
personnel  and  they  contact  us  whenever  they 
encounter  problems  with  the  code  that  they  cannot 
solve  themselves.  Our  experts  are  regularly  (2-3 


times  a  year)  visiting  CCMC  to  help  to  install  new 
code  versions  and  solve  any  problems  that  may 
arise.  As  a  result  of  this  close  collaboration, 
CCMC  is  running  SWMF  in  real  time  and  publish 
real-time  plots  on  the  web.  This  is  a  great  service 
to  the  community  and  a  very  good  test  of  the 
code.  In  fact,  this  is  a  necessary  intermediate  step 
in  transitioning  SWMF  to  Air  Force  and  NOAA 
SEC  operations. 

In  our  opinion  the  collaboration  with  CCMC  is 
extremely  beneficial  to  both  sides  and  wc  consider 
it  a  model  for  other  collaborations  (with  the  Air 
Force  and  with  NOAA  SEC). 

2.2  Selected  Science  Highlights 

As  it  can  be  seen  from  the  attached  lists  of  peer 
reviewed  publications,  invited  conference 
presentations  and  contributed  talks,  our  MURI 
project  was  a  scientific  success.  The  COSTEM 
team  published  over  100  refereed  papers  in 
leading  journals  and  monographs.  Wc  gave  over 
100  invited  presentations  at  national  and 
international  conferences  and  workshops.  Last  but 
not  least,  the  team  gave  174  contributed 
presentations  at  various  scientific  meetings. 

Instead  of  listing  all  scientific  results  (that 
would  make  this  report  quite  long)  wc  simply 
refer  to  the  publications  and  enclose  a  PDF  copy 
of  all  papers  listed  in  the  publication  list.  In  this 
section  wc  just  list  a  few  scientific  highlights  to 
illustrate  our  accomplishments.  This  list  is 
incomplete  and  somewhat  subjective,  but  it  gives 
a  cross-section  of  our  scientific  accomplishments. 

2.2. 1  Flare  Heating  [68] 

We  have  carried  out  a  theoretical  analysis  of 
the  thermal  radiation  emitted  by  large,  eruptive 
Hares.  The  analysis  is  based  on  the  configuration 
shown  in  Figure  8  which  consists  of  an  upward 
moving  magnetic  fiux  rope  with  a  vertical  current 
sheet  below.  Reconnection  at  the  current  sheet 
converts  the  magnetic  energy  of  the  plasma 
flowing  into  the  sheet  into  kinetic  energy  and 
heat.  The  analysis  assumes  that  at  least  half  of  the 
Poynting  flux  into  the  sheet  is  channeled  along 
field  lines  to  the  chromosphere  where  it  drives  an 
upflow  of  dense  plasma.  This  process  is  known  as 
a  chromospheric  evaporation,  and  it  leads  to  the 
formation  of  a  system  of  thermal  flare  loops  as 
shown  in  Figure  8, 
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The  temperatures  and  densities  resulting  from 
chromospheric  evaporation  were  calculated  using 
the  simple  evaporative  cooling  model  of  Cargill 
et.  ah  These  values  were  subsequently  used  to 
determine  theoretical  flare  light  curves  for  the 
Transition  Region  and  Coronal  Explorer 
(TRACE),  the  Soft  X-ray  Telescope  (SXT)  on  the 
Yohkoh  satellite  and  the  Geostationary 
Operational  Environmental  Satellite  (GOES),  The 
correlation  between  the  speed  of  material  ejected 
as  a  coronal  mass  ejection  (CME)  and  any 
associated  flare  is  not  straightforward*  For 
example,  it  is  possible  to  have  two  CMEs  with 
nearly  the  same  trajectories  and  speeds  but  for 
which  there  is  a  tenfold  difference  m  the  peak 
intensities  of  their  light  curves. 

The  magnetic  configuration  used  for  the 
calculation  is  based  on  a  loss  of  global,  ideal- 
MHD  equilibrium  in  a  flux  rope  that  is  suspended 
in  the  corona  by  a  balance  between  magnetic 
tension  and  compression.  Equilibrium  is  lost 
when  the  magnetic  boundary  condition  at  the 
photosphere  is  slowly  evolved  to  a  critical  point 
where  a  balance  between  compression  and  tension 
is  no  longer  possible.  When  this  point  is  reached, 
the  flux  rope  erupts  outwards  to  form  a  vertical 
current  sheet  as  shown  in  Figure  8. 

From  our  analysis  we  have  found  that  the 
fraction  of  the  released  magnetic  energy  that  goes 
into  thermal  energy  depends  strongly  on  the 
reconnection  rate  at  the  current  sheet.  As  a 
measure  of  the  reconnection  rate  we  use  Ma.  the 
AlfVen  Mach  number  at  the  midpoint  of  the  edge 
of  the  sheet.  For  MA  near  unity  only  about  15% 


of  the  total  magnetic  energy  in  the  configuration  is 
converted  into  the  thermal  energy  of  the  flare. 
The  remaining  85%  is  channeled  into  the  kinetic 
energy  of  the  ejected  mass.  However,  as  the 
reconnection  is  made  more  difficult  (e.g.  by 
reducing  the  electrical  resistivity  of  the  plasma), 
the  percentage  of  the  total  magnetic  energy  which 
is  thermalized  increases.  For  MA=  0,01,  about 
80%  of  the  energy  is  thermal  and  only  about  20% 
is  transferred  to  the  kinetic  energy  of  the  ejecta. 
These  variations  reflect  the  fact  that  as  the 
reconnection  rate  decreases  the  current  sheet 
becomes  longer,  so  that  even  though  the  rate  at 
which  magnetic  energy  flows  into  the  sheet 
decreases,  the  net  Poynting  flux  into  the  current 
sheet  increases.  As  MA  tends  to  zero  the  thermal 
energy  eventually  goes  to  zero  in  the  absence  of 
any  reconnection, 

2. 2, 2  ICME  Structure  [ 28 \  42] 

In  an  effort  to  understand  the  structure  and 
time  evolution  of  such  CME-driven  shocks  and 
their  relevance  to  particle  acceleration,  we 
investigate  the  interaction  of  a  fast  CME  with  the 
ambient  solar  wind  by  means  of  a  three- 
dimensional  numerical  ideal  MHD  model.  Our 
global  steady-state  coronal  model  possesses  high- 
latitude  coronal  holes  and  a  helmet  streamer 
structure  with  a  current  sheet  near  the  equator, 
reminiscent  of  near  solar  minimum  conditions. 
Fast  and  slow  speed  solar  wind  flow  at  high  and 
low  latitude  respectively  and  the  Archimedian 
spiral  geometry  of  the  interplanetary  magnetic 
field  is  reproduced  by  solar  rotation.  Within  this 
model  system,  we  drive  a  CME  to  erupt  by  the 
introduction  of  a  Gibson-Low  magnetic  flux  rope 
that  is  embedded  in  the  helmet  streamer  in  an 
initial  state  of  force  imbalance.  The  flux  rope 
rapidly  expands  and  is  ejected  from  the  corona 
with  maximum  speeds  in  excess  of  1000  km/s 
driving  a  fast-mode  shock  from  the  inner  corona 
to  a  distance  of  1  AU, 

We  find  that  the  ambient  solar  wind  structure 
strongly  affects  the  evolution  of  the  CME-driven 
shocks  causing  deviations  of  the  of  the  fast-mode 
shocks  from  their  expected  global  configuration. 
These  deflections  lead  to  substantial  compressions 
of  the  plasma  and  magnetic  field  in  their 
associated  sheath  region.  The  sudden  post-shock 
increase  in  magnetic  field  strength  on  low  latitude 
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Held  lines  is  found  to  be  effective  for  accelerating 
particles  to  the  GeV  range. 


This  same  simulation  demonstrates  an 
alternative  explanation  for  forward- reverse  shock 
pairs  observed  by  Ulysses  to  bound  high-latitude 
CMEs,  It  has  been  suggested  that  these  forward- 
reverse  shock  pairs  form  as  a  result  of  coronal 
mass  ejections  into  the  ambient  solar  wind,  so 
called  “over-expansion,"  Within  our  CME  model, 
wc  find  that  when  the  CME  is  greater  than  40  Rs 
from  the  Sun,  a  reverse  shock  forms  poleward  of 
the  CME  as  a  result  of  the  interaction  of  the  CME 
with  the  solar  wind.  In  front  of  the  CME,  the 
slow  wind  is  deflected  to  higher  latitude  while 
behind  the  CME,  fast  wind  is  deflected  to  low 
latitude.  These  deflected  streams  collide  to  form  a 
reverse  shock.  The  shock  pair  formed  in  this  way 
naturally  occurs  at  high  latitude  in  the  fast  wind  as 


observed,  and  explains  many  features  associated 
with  the  shock  pair. 

2.2.3  Interacting  CMEs  [2] 

Large  active  regions  often  produce  CMEs  in 
rapid  succession.  When  a  later  CME  is  faster 
than  the  earlier  one(s)  complicated  interactions 
can  take  place  in  the  interplanetary  medium 
that  fundamentally  change  the  plasma 
properties  in  the  heliosphere.  Wc  published 
several  papers  investigating  this  phenomenon 
[2,33]. 

Here  we  show  MHD  simulation  results  for 
coronal  mass  ejections  originating  from  NOAA 
active  region  9236  on  November  24,  2000. 
These  three  ejections,  with  velocities  around 
1200  km/s  and  associated  with  X-class  flares, 
erupted  from  the  Sun  in  a  period  of  about  16.5 
hr.  In  our  simulation,  the  coronal  magnetic  field 
is  reconstructed  from  MDI  magnetograin  data, 
the  steady-state  solar  wind  is  based  on  a  varying 
polytropic  index  model,  and  the  ejections  are 
initiated  using  out-of-equilibrium  semi- 
cylindrical  flux  ropes  with  a  size  smaller  than 
the  active  region.  The  simulations  were  carried 
out  with  the  SWMF. 

The  use  of  an  out-of-equilibrium  flux-rope 
model  to  initiate  the  eruptions  makes  it  possible 
to  reproduce  qualitatively  and  quantitatively  the 
features  of  LASCO  observations.  We  find  that  a 
CME  originating  from  the  Sun  center  docs  not 
necessarily  appear  as  a  symmetric  halo  in  line- 
of-sight  images.  The  initial  deflection  of  the 
ejection  is  determined  by  the  interaction  of  the 
ejection  with  the  solar  magnetic  field.  In  this 
case,  the  ejection's  speed  in  the  plane  of  sky 
can  be  an  overestimation  of  the  real  speed  of 
the  ejection  toward  the  Earth. 

In  our  simulation,  the  transit  time  to  Earth 
of  the  complex  ejecta  resulting  from  the  three 
ejections  is  10  hr  longer  than  the  observed 
transit  time  of  54.5  hr  for  the  first  shock.  A 
previous  ejection,  the  driver  of  the  shock 
observed  by  Wind  at  05:30  UT  on  November 
26,  was  not  included  in  the  simulation.  This 
previous  CME  preconditioned  the  solar  wind, 
decreasing  the  background  density  and 
increasing  the  speed.  Therefore,  its  absence  in 
our  simulation  might  contribute  to  the 
overestimation  of  the  transit  time  of  the 
leading  shock.  In  addition,  at  Earth,  the 
simulated  solar  wind  density  is  too  large  by  a 
factor  of  2,  and  the  simulated  magnetic  Held  too 
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weak  by  a  factor  of  2,  which  would  make  the 
simulated  shocks  slower  than  the  observed  ones. 

Using  different  models  for  the  solar  wind  and 
the  magnetic  flux  ropes  than  in  our  previous 
simulation  [331,  wc  find  that  the  shock 
associated  with  the  second  ejection  remains  a 
fast-mode  MHD  shock  at  all  times.  As  it 
propagates  inside  the  first  magnetic  cloud,  the 
jumps  in  pressure  and  density  across  the  shock 
are  low,  In  part  of  the  cloud,  the  jump  in 
pressure  is  found  to  be  only  a  few  percent  above 
the  reversible  compression.  This  could  explain 
why  at  Earth  there  arc  not  always  clear 
signatures  of  the  shock  passage  inside  a  cloud 
other  than  the  shorter  duration  of  the  cloud.  As 
the  initial  density  and  thermal  pressure  inside  a 
magnetic  cloud  are  low,  the  passage  of  a  weak 
shock  does  not  necessarily  lead  to  a  clear 
signature  in  thermodynamic  quantities. 

As  the  first  two  shocks  collide,  there  is  a 
temporary  increase  in  the  density  jump,  which 
can  be  above  the  compression  limit  for  a  single 
shock.  Such  a  large  compression  should  be 
observable  by  the  Heliospheric  Imagers  onboard 
STEREO,  This  could  be  the  clearest 
observational  proof  of  the  fact  that  a  shock 
propagates  inside  a  cloud.  After  the  interaction 
of  the  first  two  shocks,  the  remaining  shock  in 
front  of  the  magnetic  clouds  is  the  faster  shock 
associated  with  the  second  cloud.  Indeed,  its 
speed  and  direction  are  consistent  with  the 
speed  and  direction  of  the  second  shock  before 
the  shocks'  interaction.  Because  a  faster, 
stronger  shock  is  now  propagating  into  the 
undisturbed  solar  wind,  the  jumps  in  pressure  and 
density  arc  greater  than  the  jumps  before  the 
shocks’  interaction.  The  third  shock  becomes  a 
compression  wave  when  it  enters  the  second 
cloud  This  is  due  to  the  large  Alfven  speed 
inside  the  preceding  cloud  and  the  slow  speed  of 
the  third  shock  relative  to  the  preceding  cloud. 

2-2.4  Heliospheric  Tomography  [22] 

We  used  3D  reconstruction  technique  to 
analyze  SMEl  data  and  provide  density  in  the 
interplanetary  medium.  These  analyses  are 
currently  being  compared  with  LASCO 
coronagraph  data  and  with  interplanetary 
scintillation  data  for  selected  time  intervals  (late 
May,  2003  and  late  October  2003).  Figure  10 
shows  the  October  28,  2003  CME  reconstruction, 
and  its  comparison  with  interplanetary 
scintillation  gdevcl  observations  obtained  during 


the  same  time  interval.  Notable  in  both 
reconstructions  is  that  the  mass  derived  from  these 
two  techniques  gives  approximately  the  same 
value  and  that  this  mass  is  in  approximately  the 
same  location  in  spite  of  the  extremely  fast  shock 
that  preceded  the  CME  and  that  has  at  this  time 
reached  the  Earth.  Past  studies  have  indicated  that 
shocked  plasma  may  contain  more  small-scale 
(-200km)  turbulence  by  as  much  as  an  order  of 
magnitude  than  other  interplanetary  regions.  If 
this  were  the  case  1  PS-derived  masses  would  be 
highly  unreliable  so  this  would  negate  the  idea 
that  the  scintillation  process  could  be  used  to 
determine  reliable  bulk  densities.  Although  more 
complete  analyses  with  other  events  are  expected 
to  refine  this  study,  the  clear  indication  for  this 
very  large  event  is  that  the  differences  between 
the  IPS  and  white  light  density  measurements  are 
minimal.  The  masses  using  the  SM.E1  data  in  the 
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figure  caption  give  an  interesting  comparison  with 
LASCO  coronagraph  observations  since  LASCO 
observed  the  upper  portion  of  the  CME  to  have 
-2*1016  g  out-of-the-sky-plane  excess  mass,  and 
the  southern  portion  to  have  ~4*10U>  g  out-of-the- 
sky-plane  excess  mass. 

2.2.5  Analytic  SEP  Model  [46] 

A  calculation  of  the  proton-excited  wave 
intensity  upstream  of  a  stationary  planar  shock  has 
been  carried  out.  The  new  calculation  relaxes  the 
assumption  made  in  the  previous  analysis  that  the 
growth  rate  of  a  given  wave  with  wave  number  * 
is  dominated  by  the  resonant  protons  with  the 
lowest  possible  energy.  Although  this  assumption 
is  reasonable  close  to  the  shock,  it  is  not  so 
reasonable  farther  upstream  because  the  upstream 
region  is  increasingly  dominated  by  higher  energy 
particles  which  can  more  easily  escape  from  the 
region  near  the  shock. 

The  calculation  reveals  that  the  wave  intensity, 
which  is  ~k*-6  at  small  *,  transitions  to  being  ~*"2 
at  large  k  with  a  transition  wave  number  *o(z) 
which  decreases  as  an  inverse  power  of  increasing 
distance  from  the  shock.  Here  p  (=  3x  ^Cx_l))  *s 
the  standard  power-faw  spectral  index  for  a 
stationary  planar  shock,  where  x  is  the  shock 
compression  ratio.  This  form  has  important 
consequences  for  the  variation  of  ion  composition 
throughout  the  event  including  the  ions  which 
arrive  promptly  at  the  observer  with  high 
streaming  anisotropy.  The  low-*-  power  law, 
which  resonates  with  higher  energy  ions,  leads  to 
enhancement  of  the  heavy  ions  (larger  A/Q)  with 
increasing  distance  upstream  of  the  shock* 
However,  the  large-*  power  law,  ~*~\  which 
resonates  with  lower  energy  ions,  docs  not 
fractionate  between  species  with  different  A/Q. 
The  resulting  compositional  variation  seems  to 
account  for  that  observed  in  many  events. 

We  have  also  investigated  the  origin  of  the 
large  variations  in  the  Fe/O  ratio  at  high  energies 
observed  between  different  solar  energetic  particle 
(SEP)  events.  The  origin  appears  to  arise  from  the 
magnetic  obliquity  of  the  shock.  The  ion  seed 
population  for  injection  at  the  shock  consists  of 
both  solar  wind  ions  and  ambient  energetic 
particles.  Quasi-parallel  shocks  (with  upstream 
magnetic  field  primarily  parallel  to  the  shock 
normal)  have  a  lower  energy  threshold  for  ion 


injection  than  do  quasi-perpendicular  shocks* 
Thus,  quasi-parallel  shocks  accelerate 
predominantly  solar  wind  ions,  whereas  quasi- 
perpendicular  shocks  accelerate  predominantly 
ambient  energetic  ions,  which  include  material 
from  earlier  impulsive  events  rich  in  heavy  ions. 
This  feature  of  shock  acceleration  appears  to 
account  for  the  extreme  compositional  variations 
often  observed  between  events  at  high  energies. 

2*2*6  Modeling  the  Inner  Magnetosphere 
[52] 

The  inner  edge  of  Earth's  plasma  sheet  tends  to 
shield  the  region  earthward  of  it  from  the  main 
force  of  the  magnctosphcric-convcction  electric 
field,  but  both  models  indicate  that  changes  in  the 
rate  of  convection  can  cause  temporary 
penetration  of  the  shielding.  A  southward  turning 
of  the  IMF  causes  undershielding  and  penetration 
of  an  eastward  electric  field  across  the  day  side 
ionosphere,  raising  the  F-layer,  Major  magnetic 
storms  generate  StorrrvEnhanced-Dcnsity  (SED) 
events:  GPS  observations  of  Total  Electron 
Content  show  great  plumes  of  enhanced  density 
that  sweep  across  the  afternoon  sector  of  the  mid- 
and  high-latitude  ionosphere,  leading  to 
substantial  space  weather  effects,  including 
equatorial  bubbles  and  disruption  of  the  GPS 
system*  SED  events  apparently  result  from  the 
dayside  eastward  penetration  electric  field,  w'hich 
makes  understanding  of  penetration  fields  a 
priority  issue  from  a  space-weather  point  of  view. 

The  strength  and  duration  of  the  prompt- 
penetration  electric  fields  depend  on  the  changes 
in  the  overall  magnetic  configuration  of  the 
magnetosphere  that  occur  in  response  to  changes 
in  solar-wind  driver.  Thus  proper  first-principles 
modeling  of  the  phenomenon  requires  a  code  that 
couples  both  inner  and  outer  magnetosphere  to  the 
ionosphere.  BATS-R-US/RCM  is  the  first  such 
model,  so  we  arc  applying  it  to  the  prompt- 
penetration  problem. 

The  code  was  run  for  four  hours  with  steady 
northward  IMF  (B*=5  nT,  p=5  amu/cm3,  Vx=-400 
km/s).  The  southward  turning  (to  B^-5  nT)  hit 
the  sunward  boundary  of  the  simulation  box 
(X=32  Re)  at  12:00  UT,  The  effect  reached  the 
dayside  magnetopause  about  12:18  UT,  Figure  1 1 
shows  a  sequence  of  equipotential  patterns  in  the 
northern  ionosphere.  Note  that  the  disturbance 
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effect  starts  approximately  simultaneously  in  the 
entire  ionosphere  and  gradually  builds  in  strength, 
but  there  are  interesting  subtleties.  The 
penetration  field,  as  evidenced  by  equipotentials 
that  penetrate  to  the  equator,  occur  simultaneously 
with  similar  strengths  on  both  day  and  night  sides. 
The  dayside  equatorial  electric  field  is  eastward, 
tending  to  lift  the  F-layer 


3  Lessons  learned 

In  our  opinion  this  was  a  very  successful  and 
enjoyable  project.  It  met  its  goals,  it  produced 
good  science,  it  contributed  significantly  to 
graduate  and  postdoctoral  education,  and  it  helped 
to  propel  the  careers  of  several  very  talented 
scientists. 

There  arc  several,  mainly  management  related, 
lessons  that  we  learned  in  the  course  of  this 
project.  In  this  short  section  we  will  list  some  of 
the  important  lessons  we  learned. 


One  of  the  important  lessons  is  that  the  PI  of  a 
MUR  I  consortium  has  very  limited  leverage  with 
his  team  members  located  at  other  institutions. 
Most  of  the  prestige  and  visibility  associated  with 
the  project  is  focused  at  the  PI  institution  and  the 
MURI  project  is  usually  not  the  highest  priority  of 
the  collaborating  partners.  They  tend  to  do  the 
research  what  they  would  be  doing  anyway,  using 
the  additional  MURI  funds  to  supplement  their 
existing  funding  base.  Lesson:  The  PI  { ami  the 
program  monitor)  should  not  expect  team 
members  at  other  institutions  to  do  anything  that 
is  out  of  the  ordinary  for  them. 

A  corollary  to  the  previous  lesson  is  that  when 
a  MURI  proposal  is  submitted,  the  overall  success 
or  failure  must  depend  on  the  performance  of  the 
PI  institution  and  it  is  a  mistake  to  allow  any 
partner  institutions  to  be  a  “single  point  failure" 
element. 

Graduate  students  arc  usually  very  smart 
individuals,  but  they  arc  still  at  an  early  stage  of 
their  professional  development.  They  still  are 
learning  what  is  important  and  what  can  be  done 
later  or  at  a  slower  pace.  One  of  the  important 
things  they  usually  arc  still  struggling  with  is  time 
management.  As  a  consequence  of  these  factors  it 
is  not  a  £)*od  idea  to  give  them  tasks  that  arc 
“mission  critical"  for  the  success  of  a  large  MURI 
project.  Lesson:  Graduate  students  should  be 
given  challenging  tasks ,  hut  these  tasks  need  to  be 
somewhat  detached  from  the  mainstream  of  the 
project,  and  should  not  be  time  critical 

MURI  projects,  by  their  nature,  are  large 
efforts  that  require  full-time  efforts  by  a  team  of 
well  qualified  scientists  who  arc  able  and  willing 
to  work  as  a  team.  Since  MURI  projects  primarily 
fund  university  research  groups,  it  is  important  to 
have  a  core  of  scientists  who  can  devote  most  of 
their  time  and  energy  to  the  success  of  the  project. 
In  practice  this  means  that  one  cannot  have  a 
successful  MURI  team  that  consists  only  of 
professors  and  graduate  students.  Professors  have 
far  too  many  other  responsibilities  and 
commitments  and  graduate  students  are  not  ready 
yet  to  carry  out  the  mission  critical  part  of  the 
work.  Lesson:  Successful  MURI  teams  need  to 
have  a  professional  core  team  of  postdocs  and 
research  scientists  whose  primary  job  is  the 
successful  completion  of  the  project.  Professors 
need  to  lead  the  team  and  provide  intellectual 
leadership,  and  graduate  can  participate,  hut  the 
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most  critical  tasks  must  be  performed  by 
pro fess  itma  /  sc  i  en  fists : 

MURI  projects  typically  involve 
interdisciplinary  teams  that  need  to  work  together 
efficiently.  Our  team  established  regular  weekly 
meetings  that  greatly  helped  us  to  keep  the  team 
focused  and  stay  "on  the  same  page/’  In  addition, 
we  tightly  integrated  the  team  of  professional 
scientists,  and  software  development,  algorithm 
design  and  physics  application  was  done  by  a 
single,  integrated  team.  On  the  negative  side  this 
significantly  increased  the  daily  frustration  level 
and  friction,  but  this  model  resulted  in  a  quick  and 
efficient  resolution  of  emerging  obstacles.  Lesson: 
Integrated  reams  are  very  efficient,  but 
occasionally  painful.  Nevertheless,  they  are  a 
major  factor  in  success. 

One  of  the  very  important  factors  in  our 
success  was  that  frequent  communications 
between  the  PI,  the  team  and  the  AFOSR  Program 
Director.  This  good  relationship  enabled  us  to 
understand  the  priorities  and  expectations  of  the 
funding  agency,  while  the  Program  Director  had  a 
good  understanding  of  our  progress  and  problems. 
This  enabled  the  Program  Director  to  have  the 
AFOSR  management  engaged  and  interested  in 
this  MURI  project.  Lesson:  Good  communication 
with  the  funding  agency  is  essential. 
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